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ABSTRACT 

Stars on eccentric orbits around a massive black hole (MBH) emit bursts of gravitational 
waves (GWs) at periapse. Such events may be directly resolvable in the Galactic centre. How- 
ever, if the star does not spiral in, the emitted GWs are not resolvable for extra-galactic MBHs, 
but constitute a source of background noise. We estimate the power spectrum of this extreme 
mass ratio burst background (EMBB) and compare it to the anticipated instrumental noise of 
the Laser Interferometer Space Antenna (LISA). To this end, we model the regions close to 
a MBH, accounting for mass-segregation, and for processes that limit the presence of stars 
close to the MBH, such as GW inspiral and hydrodynamical collisions between stars. We find 
that the EMBB is dominated by GW bursts from stellar mass black holes, and the magnitude 
of the noise spectrum is at least a factor ^10 smaller than the instrumental noise. 

As an additional result of our analysis, we show that LISA is unUkely to detect relativistic 
bursts in the Galactic centre. 
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1 INTRODUCTION 

A new opportunity to study stellar processes near massive black 
holes (MBHs) arises with the anticipated detection of gravi- 
tational waves (GWs) by the Laser Interferometer Space An- 
tenna (LISA). LISA will be a space-based detector in orbit 
around the Sun, consisting of three satellites five million kilo- 
metres apart. It will be sensitive to GWs in a frequency range 
10"'* Hz < / < 10"- Hz. An important source of GWs for LISA 
is the in spiral of compact ob j ects onto MBHs in galactic nu- 



clei (e.g. Illils & Bendeilll995l: ISigurdsson & Reeslll997| ; 
2OO2I: lFreitag|l2003l; iHopman & Alexandeil I200I l2006a 



Amaro-Seoane et '^ 120071 : see lHopmanll2006l for a review). These 



are sources on highly eccentric orbits with periapses slightly larger 
than the Schwarzschild radius rs = 2GM,/c^ of the MBH, where 
M, is the mass of the MBH. The star dissipates energy due to 
GW emission, and as a result spirals in. Such extreme mass ra- 
tio inspirals (EMRIs) can be observed by LISA to cosmological 
dist ances if the orbital period of the star is shorter than P < 10^ 
sec jpinn & Thorne' '200(]| : iBarack & Cutie3l2004 iGair et aljr2004 
iGlampedaidsi ,2005i) . LISA will detect hundreds to th ousands of 
such captures ov er its projected 3-5 yr mission life time l lGair et al] 
l2004lGaiil200g) . 

For most of the inspiral the emitted GWs are not observable by 
LISA. These GWs give rise to a source of confusion noise, possibly 
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obscuring other types of GW sources. The shape a nd overall magni- 
tude o f this EMRI background has been studied bv lBarack & Cutleii 
( I2OO4I) . They do not study the dynamical requirem ents of inspiral , 
but sc ale their result with possible EMRI rates of lFreita3 l l200ll . 
l2003h . Therefore noise of stars that do not eventually spiral in is 
not included. 

iRubbo et al. I ( I2OO6I) show that stars on long periods of a 
few years and nearly radial orbits that carry them near the 
Schwarzschild radius of the MBH, emit bursts of GWs that for the 
Galactic centre will give a signal to noise ratio larger than 5. Taking 
into account processes determining the inne r radius of the density 
profile as well as mass segregation effects, iHopman et alj l l2007h 
find that stellar mass black holes (BHs) have a burst rate of order 
1 yr"', while the rate is < 0.1 yr"' for main sequence stars (MSs) 
and white dwarfs (WDs). 

Individual bursts from star-MBH fly-bys may be detected 
from our own Galactic centre, but it is unlikely they will be ob- 
served from other galactic nuclei. However, the accumulation of all 
bursts of all galaxies in the universe gives rise to an extreme mass 
ratio burst background (EMBB). Since the energy emitted per event 
is much higher for inspirals compared to fly -bys, but the event rate 
is much lower ( [Alexander & Hopman|[2003l) . it is not a priori clear 
which of these dominates the confuse background. In this paper we 
study the contribution of these fly-bys to the GW background. In ^ 
we derive an analytical expression for the event rate of GW bursts 
in a galactic nucleus. This model diverges at several boundaries, 
and we consider in 92.21 what the physical processes are that de- 
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termine the range within which our model i s valid. Our formalism 
is similar to that of the galactic burst rate of lHopman et alj l l2007h . 
We derive an expression for the emitted energy spectrum of GW 
bursts in a single galactic nucleus in ^ A relation for the EMBB 
from the accumulation of GW bursts from all redshifts observed 
by LISA is derived in 33.31 In j4]w e compare the resulting EMBB 
to the instrumental noise of LISA ( lLarsonll200lh . and to other as- 
trophysical backgrounds. We conclude in ^and discuss what our 
model implies for the possibility of directly detecting relativistic 
bursts of GWs from our Galactic centre. 



2 THE GRAVITATIONAL WAVE BURST RATE AND 
SPECTRUM FROM A SINGLE GALACTIC NUCLEUS 

In this section, we present a simple model for a galactic nucleus 
and derive the burst rate for that model. We discard parts of phase- 
space in our model where stars are unlikely to exist due to various 
processes. 



2.1 The gravitational wave burst rate 

A MBH dominates the potential inside the radius of influence 
GM. 



Ill ■- 



(1) 



where a denotes the stellar velocity dispersion far away from the 
MBH. The MBH mass is empirically related to the velocity disper- 
sion of the bulge by 



M. 



(2) 



,200 km s" 

(e.g. lMerritt & Ferraresel200ll ; lTremaine et alj2002h . Inside the ra- 
dius of influence, stellar orbits are assumed to be Keplerian. 

Galactic nuclei with M. < 10^ M ;, have relaxation times less 
than a Hubble time and are dynamically relaxed, while for very 
massive MBHs, the distribution will depend on the initial con- 
ditions. However even in this case the inner region of the cusp 
may be of power-law form due to the effect caus ed by a slow 
growth of a black ho le inside a stellar system (see e.g. lYoun jl980l ; 
IOuinlanetal]|l995h . We will for simplicity assume that our mod- 
els also apply for M, > lO^M;,. The contribution of such MBHs is 
only marginal for the EMBB, and this assumption does not affect 
our result. 

We assume the distribution within r/, to be spherically sym- 
metric in space and approximately isotropic in velocity space. 
The radial density profile is approximated by a power law, v = 
r-"" (e.g. IPeeblesIl 19721 : iBahcall & Woljl 19761) . where the slope 
OtM depends on the stellar mass M^,^, due to mass segregation 
ISDitzej ( fl987l) , larger masses have larger 0!^/. The number of stars 
nM{a,fi) da dj~ of type M in an elemsnl {a, a + da),{J-, fi + 
dfi) is then given by 



nM(a,j'') da dj^ = (3 — CCm) 



CmN,, 



JHa\ 



dadJ-, (3) 



where we assume A'/, = M, /Mq, . The total number of stars of type 
M within r/, is defined as CmN^- 

We calibrate the parameters Cm using solutions to 
Fokker-Planck equations a pplied to the Galactic centre 
( [Alexander & Hopmad l2009h . such that the number of stars 
within 0.01 pc agrees with those simulations. Since the density 



distribution of a species M is not a strict power law in the Fokker- 
Planck models, this implies that our density profile is a good 
approximation only at small radii. Since the burst rate is dominated 
at those small radii (see equation (7)), the small mismatch at larger 
radii will not strongly affect our conclusions. 

A useful parametrisation of an orbit is the frequency a star 
would have on a circular orbit with radius of its periapse rp. 



CO,, 



GM. 



(4) 



Note that for eccentric orbits cOp ^ l/P where P{a) = 
2K{a^ /GM,yl'^ is the orbital period. The number of stars 
n{a, COp) da dcOp in an element (a, a + da),{cOp, cOp+dcOp) is then 



nM[a,co 



■pi 



:(3-«m)^(GM.)V3 

3rf 



l-au 



CO, 



-5/3 



(5) 



where it was assumed that e ~ 1. 

A star on a highly eccentric orbit emits almost all of its GW 
energy at periapse in a single GW burst. The burst rate per loga- 
rithmic unit of the semi-major axis and frequency at which stars of 
species M have a bursting interaction with the MBH is given by 



d^FM 



dlna dlnco,, 



^^(3-«m)^(GM.) 

P{a) 3rh 

:5Jx 10"^ yr"' {3-aM)CMM 



)2—aM 



(6) 



0!„/2+l/3 - l/2-a„ -2/3^ 



where we convert to dimensionless units defined as M, = M,/Mo; 

a = ajrh„ \ cOp = cOp / cOo , and Mo = IO^Mq,; r/,,o = GMo/CTo ^ Ipc; 
fflo =cV(GMo) Ri0.2Hz. 

2.2 The inner region of the stellar cusp 

Since we describe the profile of the galactic nucleus by power 
laws, integration over orbital space leads to formal divergences, 
and it is therefore important to determine the boundaries of va- 
lidity of our model. In the next section we consider a num- 
ber of processes that can determine the inner edge of the stellar 
cusp. The burst rate is approximately proportional to dTM/dlna o= 
^jl-au (-ggg equation (7)). As typical values for a are in the range 
1.4 — 2 (see e.g. |Hopman & Alex ander 2006a; Freitag et al. 2006|; 
[Alexander & HopmanI floogi) . the GW energy spectrum is domi- 
nated by GW bursts from stellar orbits at small radii and therefore 
we focus on the stellar dynamics in the vicinity of a MBH. 

2.2.1 Gravitational wave iiispiral 

As a star loses energy repeatedly with every periapse passage due to 
GW emission, it spirals in. Unresolv able EMRIs cons t itute a back- 
ground noise for LISA studied by iBarack & CiitleJ | |2004) . Our 
study instead focuses on the subset of extreme mass ratio events 
that do not spiral in, and GW emission from EMRIs must be ex- 
cluded. We note that the dynamics of the progenitors of EMRI noise 
and burst noise is very different. The part of phase space where 
stars spiral in is at any particular time typically not populated, so 
that, e.g., equation (Q does not apply. 

The amount of energy A Eqw that is lost in one period is given 
bv lPeters & Mathewsl ( 1 19631) as 
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AE, 



GW 



7= /(e) 

5V2 M. V 



-7/2 



1 X 10^9 ergM2^.'/3 -7/3^ 



(7) 



where for highly eccentric orbits /(e) ~ 0.39. We define the char- 
acteristic time for inspiral f,- as the time it takes the initial orbital 
energy Eq to grow formally to infinity, 



ti 



Ie„ de/dt Jen 



de 



A£Gw//'(fl) 



2.9 X 10^ yrM^' M. 



(8) 



While the orbit decays, two body scattering by other stars per- 
turbs it, changing the orbital angular momentum by order of itself 
in a timescale tj. The magnitude and direction of the step in angu- 
lar momentum is a random walk process, and typically in a relax- 
ation time the step is of the size of the circular angular momentum 
Jc = {GM.aYl'^, and hence 



tj- 



^rela 



In this expression r^elax is the relaxation time, 



^relax ^relax 



1/2 



(9) 



(10) 



where Qeiax is a numerical constant that also absorbs the Coulomb 
logarithm, o is the velocity dispersion far away from the MBH, and 
nil °= M, ~ is the density of stars at th e radius of influence. 

Due to strong mass-segregation l i Alexander & Hopmanll200^ ; 
iKeshet et al at small radii the relaxation time is dominated 

by BH interactions, while at large distances the MSs dominate the 
relaxation rate. The true relaxation time is a combination of the 
BH and MS relaxation time, and in order to keep the calculation 
more tractable, we approximate it by a power law. As the burst rate 
is critically sensitive to the inner cut-off of the density profile, we 
chose the constant Qeiax to calibrat e the relaxation time at a sma ll 
radius of 0.01 pc for the model in [Alexander & HopmaiJ ( l2009h . 
For those values, we find Qeiax = 0.02. 

With this scaling the relaxation time and scatter time become 



W = 2.7 Gyr^M.«l/2; 

^ 0.02 ' 



(11) 



\i tj > ti, stars spiral in faster than they are replenished by 
other stars. In the corresponding region of orbit space, any bursting 
star is quickly accreted, so this region is typically empty (see the 
appendix for further details). The fact that this part of phase-space 
is not populated has also consequences for the possibility of de- 
tecting relativistic bursts in the Galactic centre; see ^for further 
details. Solving fy = f; for (Hp gives an inner cut-off of 



a, ^ 1.1 (^)-3/5 ^ -3/5 ^ -13/10 3/5^ 
^ 0.02 ^ 



(12) 



We assume that stars with (Bp > ffl; always spiral in and we 
neglect these stars (n oise for these stars was accounted for by 
iBarack & Cu"tie3l2004t) . while stars with (bp < (B/ never spiral in, 
and we take these stars into account as bursting sources. 



2.2.2 The direct capture loss-cone 

All stars with an angular momentum less than the last stable orbit 
■^LSO = 4GM, I c (the "loss-cone") are pulled into the MBH. Since a 
star in this region of orbit space is removed in a dynamical time, we 
do not consider stars in this region. Therefore stars only contribute 
GWs to the EMBB if 



(bp < a),c = 4.41 X 10"^ M.\ 



(13) 



In reality the empty region inside the loss-cone will 
also affect the regions n earby in angular momentum space 



jLightman & Shapiro! 1 1977 ), although th i s is p a rtially erased by 
resonant relaxation dRauch & Tremaind 1 19961 : iRauch & Ingalil 
Il998l) . Here we assume that orbits outside the loss-cone are fully 
populated. 



2.2.3 The tidal disruption loss-cone 

Tidal effects can play a role when the star is not compact. The 
tidal forces from the MBH on the MS disrupt those MSs whose 
orbits carry them within the tidal radius r, ~ (2M,/M^)'/-'i?^ of 
the MBH, where R^, denotes the radius of the MS. This effect leads 
to a depletion of MSs when r < rt, called the loss-cone for tidal 
disruptions. We only consider MSs if 



(B„ < (B, = 2.18 X 10 



_3 iMA'^' (R.^-"' 



Mr., 



R. 



(14) 



Note that the break up frequency a)t depends solely on the charac- 
teristics of the MS. 



2.2.4 Hydrodynamical collisions 

Close to a MBH, MSs are likely to undergo mu ltiple collisions in 
their lifetime (see e.g. ISpitzer & Saslawi[l96^ . The rate 

Tcoll at 

which a MS with radius i?^ has grazing collisions can be estimated 
by 



3-0!ms Nh_ 
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CM. 



(15) 



Multiple collisio ns can lead to the disruption of the MS. 
iFreitag & Ben3 l l2005l) showed that on average a MS is disrupted 
after A^coU ~ 30 collisions. This implies that MSs are saved from 
disruption by collisions within a Hubble time, if their distance to 
the MBH is larger than 



a > flcou = 1.6 X 10 



30 

A'coll 



10/19 



20/19 



Mj , 



(16) 



where it was assumed that Kms = 1.4 as in our main model. MSs 
are contributing GWs to the EMBB if their orbits obey d > flcoU- 

2.3 Model of a galactic nucleus 

To summarise, the method to compute the GW burst rate of a single 
galactic nucleus is as follows. We consider four distinct species of 
stars, that is IOMq BHs, 0.6Mq WDs, 1.4Mq neutron stars (NSs), 
and \Mq MSs. Their density profile i s given by a slope of asH = 2 , 
«WD = 1-4, Kns = 1-5, Kms = 1.4 { [Alexander & Hopmarill2009h . 
We assume that the enclosed number of MS stars within r;, is 
Nj, — M,/M-.,. We calibrate our model such that for the Galactic 
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ce ntre, the number of st a rs wit hin 0.01 pc is equal to that found 
in [Alexander & HopmanI ( |2009|) , to find Cbh = 8 x 10"^, Cwd = 
0.09, Cns = 0.01, Cms = 1- From their models, we find that the 
relaxation time at 0.01 pc in the Galactic centre is 0.78 Gyr, leading 
to a calibration factor of C^e\m — 

0.02. 

We also consider a model suggestecPlbv lO'Learv et al" U2OO8I) . 
In this model, the density of the BHs at the radius of influence 
is higher by a factor 10 at r/, compared to the main model, and 
the masses of the BHs are higher. Such models may be rele- 
vant for Galactic nuclei because of empirical hints that the ini- 
tial mas s function is much flatter than usual the Galac tic centre 
(see e.g. lNavakshin & Sunvaevl2005l:lManess et al.l2007l) . We note 
that there are no dynamical con straints on the amount or mass of 
BHs in the Galactic centre (e.g. iMouawad et al]|2005l : iGhez et al] 
l2008l) . Since in systems with many, more massive BHs mass segre- 
gation is less pronounced, the scaling for the density profile is not 
straightforward; for example, the number of BHs enclosed within 
0.01 pc is only a factor ~ 2 highe r than in our main m odel. Us- 
ing the relaxation times given in lO'Learv et J] ( |2008|) . we find 
that for this model Mbh = ISM©; Cm = 1-6 x 10^^; Kbh = 2; and 
Qeiax = 3x 10-3. 

When using these models, we exclude certain orbits due to 
the four processes described in the previous section. Stars are only 
contributing GWs if the inequalities l ll2t . J13t . l ll4t and l ll6t are 
satisfied. As an illustration we show the orbit space around a MBH 
of M. = 10*' Mq for a BH of = 10 Mq in figure ([TJi) and for a 
MS of M-t: = 1 Mq in figure (T]?). The boundaries of the populated 
parts of the orbit space of compact objects are caused by two pro- 
cesses; the inspiral process and the loss-cone. The cut-offs through 
the orbit space of MSs are caused by all four processes of 92.21 We 
note that the processes that limit the GW burst rate depend on the 
MBH mass. 



3 GRAVITATIONAL WAVE BACKGROUND 

3.1 Gravitational wave energy in a single burst 

For an eccentric Newtonian orbit, the amount of energy emitted in 
GWs in one period AE^w [see Eq. Q] and the distribution of this 



to be constant throughout the history of the Universe 

^7 



dn. 10' ^ _3 

= Gpc 

dM. M. ^ 



10^ Mr., < M. < 10*^ Me 



(18) 



The upper boundary comes from a sharp drop in the observed MBH 
number density around M. 10^ Mq. Due to a lack of information 
on MBHs of M, < 10^ Mp,, the lower boundary is arbitrarily cho- 
sen to be 10^ Mi7,. Although at frequencies between 10^^ and 10^^ 
Hz, MBHs of the lowest masses contribute most to the EMBB, de- 
creasing the lower bound only mildly affects the EMBB. The en- 
ergy density of GW bursts emitted by stars of species M is 

''^^'^ AEGwS{f-2c0p)^dM.d(0i,da. (19) 



da dcOn 



dM. 



3.3 Observed gravitational wave energy density 

The GW spectrum observed by LISA is the accumulation of the 
GW radiation of all galaxies in the universe. The emitted spectrum 
<S'{fem)dfem is the rate per unit proper time and per unit co-moving 
volume at which GW energy in the frequency range (/em' /em + 
rf/em) are emitted. The total GW energy density (dp/df)df today 
with frequencies / = /em/(l +z) in the range (f,f + df) is then 



dp 

df 



df 



dfem / dz 



dt S{U 



"dz 1+Z 



f dt , 
I dz-S[f{\+z)]. 



(20) 



For < z < 2, we follow iBarack & CutleJ ( |2004|) in approximat- 
ing f ~ /() (1 + z)-' '^, where a spatially flat Friedmann-Lemaitre- 
Robertson- Walker Universe is assumed with = 0.70 and Q.m = 
0.30 and the Universe's current age is Iq = Qi.96AHq^ = 1.39 x 
10'*Vi7() yr For sources in the range < z < 2, this is accurate 
to within ~3 per cent. The approximation is justified since the 
EMBB is dominated by the nearby universe and contributions from 
1 < z < 2 accounts for just 5-8 per cent of that of < z < 2. 



energy over frequencies dQ/d f is studied by IPeters & Mathewsl 3.4 Gravitational wave background for LISA 



( 119630 . The emitted power peaks at a frequency / that is twice 
the circular frequency at periapse cOp. To simplify our analysis we 
assume that all energy ABgw is emitted at / = 



: 2a)p, therefore 



dQ 

d.f 



(17) 



The energy density of an isotropic background of individually un- 
resolvable GW sources is related to the spectral de nsity in the LISA 
detec tor by 5gw(/) = AG/{Kc^f-) dpu/df l lBarack & CutleJ 
l2004h . or 



In reality, the spectrum is quite broad, but the exact shape of the 
spectrum is not important, since we sum over a very large number 
of bursts. 



3.2 Gravitational wave energy density 

We next integrate t his spectrum over M, with the space number 
density of MBHs of lAllen & Richstond ( |2002|) . which we assume 



•Sgw(/) = — ^ 



J J.* J J dt dn, d Fm 

„ ^ , dz dM, dcOp da 

Kc-^f- J ' dzdM,dadO)n 



AEcwS[f{l+z)-2coj 



(21) 



where equations ( 1191 ) and l l20t were used. The integral is taken over 
y representing that part of orbit space that is populated according 
to SO 



' lO'Learv et alj ilOOSi) consider a wide range of models, including models 
wit h steeper mass functio ns. We will refer to the model considered here as 
the iO'Learv et~ai] ^20081) model" for brevity. loi^eai-y et alJilOOSi) assume 
a range of BH masses, which we collect here into a single averaged BH 
mass. 



4 RESULTS 

For our main model, the resulting EMBB (/5gw)'^'^ of the four 
species of stars is shown in figure ([2j. The EMBB is dominated by 
contributions from BHs. This is mainly due to two reasons. First, a 
single burst emitted by a BH is more energetic than a burst emitted 
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Figure 1. Boundaries in orbit space around a MBH of M, = IO^Mq for a BH of M^, = lOM., in Fig.[TJi and of a MS of M^, = IM.-.j in Fig.[T]). Stellar orbits 
are populated in the shaded area only. The figures are representative of all relevant stars and MBHs where the cuts scale according to Eq. lilt . U3K )14t and 
il6\ . The red line indicates the boundary in orbit space resulting from the condition rp ^ a. The result is not very sensitive to the inner cutoff; imposing a hard 
cutoff at 10"^ pc in figure (a) decreases (/Sqw)'^^ by ~ 50 per cent. 
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Figure 2. Co mparison between the EMBB, instru mental noise of LISA 
)Larsonll200ll) , the galactic WD binary background l Larsonl l200ll) and the 
BH EMRI (inspiral) background iBarack & Cutler 2004)) with an inspi- 
ral rate of 10^'yr^'. The EMBB is smaller than the other types of GW 
noise. The EMBB due to only BHs, WDs, NSs or MSs is shown in sym- 
bols and contribution s from BHs dominate the spectrum. The EMBB of the 
lO'Learv et all )2008l) model is also shown. Even though it is significantly 
higher than the EMBB in our preferred model, it is still lower than the in- 
strumental noise of LISA. 



The EMRI background is caused by BHs that spiral in to the 
MBH. These sources are excluded in our calculations by eliminat- 
ing inspiral orbits, see 92.2.11 The expression for the EMRI back- 
ground was taken fromlBarack & Cutleil |2004|) and scales with the 
inspiral rate. From lHopman & Alexandeil ( 2006bl) we take an inspi- 
ral rate of 1 x 10^^ yr^' leading to a EMRI background at a level 
comparable to the instrumental noise of LISA. 

In our favoured model the EMBB is well below the instrumen- 
tal noise of LISA. At closest approach the EMBB (/5gw)'''~ is a 
factor ~ 10 below the instrumental noise of LISA, and a factor ^ 
19 lower than the confusion noise of Galactic WD binaries. The 
EMBB is also lower than the EMRI background in the interesting 
frequency range where the noise from EMRIs is stronger than that 
of the in strument itself. 

The lO'Learv et al. liooi) model is also displayed in figure 
In this model there are more, and more massive BHs, which raises 
the EMBB due to an increase in the number of BH bursts and an 
increase in the relaxation rate protecting stars from inspiral. Fig- 
ure 12) shows the EMBB is raised by a factor 6 on average with 
respect to our main model. However, even with this stellar distribu- 
tion the EMBB is smaller than the instrumental noise of LISA by a 
factor ^ 1 .5, and a factor ~ 3 including the noise of galactic WD bi- 
naries. Overall the sensitivity of the EMBB on the individual model 
parameters is not strong enough to raise the EMBB above the in- 
strumental noise of LISA, without increasing Mbh by a factor 10 
and other parameters by even more. 



5 SUMMARY AND DISCUSSION 



by the other species of stars in a similar orbit. Second, the BH dis- 
tribution is steeper due to mass-segregation which leads to a higher 
burst rate, see Eq. Q. Note that the EMBB due to MSs is cut off 
at ~ lO^-* Hz, because tidal effects prohibit the existence of non- 
compact objects at small distances from the MBH. 



Stars whose orbits carry them near the Schwarzschild radius of a 
MBH emit gravitational waves in the LISA frequency band. If the 
orbit of the star is very eccentric, the star emits bursts of GWs each 
peri-centre passage without necessarily spiraling in. Most extra- 
galactic bursts are not individually resolvable and hence will con- 
stitute a gravitational wave background. 
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For our favoured model the EMBB is dominated by GW 
bursts from BHs, and is a factor ~ 10 lower than the instrumental 
noise of LISA. Including the WD Galactic background increases 
the insignificance of the EMBB to a factor ~ 19. Even for the 
lO'Learv et al.l ( |2008|) model, which has a much flatter mass func- 
tion, the EMBB is smaller than the instrumental noise of LISA by 
a factor ~ 1.5 (or a factor ~ 3 including the WD Galactic back- 
;round). The EMBB i s also well below the EMRI background 
Barack & Cudedl2004h . We conclude that the detection of GWs 
by LISA will not be hindered by a background of bursting sources. 

As an additional application of our model, we revisit the possi- 
bility of direct detection o f GW bursts i n the Galactic centre, due to 
stars on relativistic orbits. lYtines et alj ( |2008|) have investigated the 
impact of including relativistic corrections to the description of the 
star's trajectory. The degree to which the relativistic corrections are 
important depen ds on the star's o rbit and is largest for small peri- 
centre distances. lYunes et al.l ^200 8) find that orbits with peri-centre 
velocities \vp\ > 0.25c account for app roximately half of th e events 



within the orbit space considered b v iRubbo et al. (120061) around 
a MBH of 3.7 x lO*" Mq. However. iRubbo et alj i2006h incorpo- 
rate tight, possibly relativistic orbits that are unpopulated accord- 
ing to our assumptions of the boundaries in orbitspace of g2.2l (see 
also the appendix). Assuming our main model o f a galactic nu- 
cleus, a similar MBH mass as lRubbo et al. U2006h and integrating 
over all frequencies that can be detected by LISA (See Eq. (4) in 
iHopman et al.l | |2007|) with a signal to noise of 5), relativistic orbits 
account for ~ I per cent of all events within the orbit space con- 
sidered by our main model. We conclude that LISA is unlikely to 
detect relativistic bursts from the Galactic centre. 



The form of this co efficient is such that it reproduces the 
iBahcall & Wolj i ll 9761) solution with distr ibution function /(x) o= 
x'/^; the prefactor is for consistency with iHopmanI ( l2009h . (Here 
positive energy implies that the star is bound to the MBH). 

Angular momentum diffusion — Let the angular momentum j 
be scaled by the circular angular momentum, such that by definition 
< j < 1. Let / perform a rand om walk such that the diffusion 
coefficient is (see lHopmaril ( l2009l) for the prefactor) 



Di 



:6.6x'/4. 



(A2) 



Boundary conditions — Let the actual angular momentum of 
the loss-cone be independent of energy, such that if it is scaled by 
the circular angular momentum, it becomes 



jlc{x) ■■ 



{-) 



1/2 



(A3) 



with Xmax some maximal energy. 

Energy loss by GWs — Let the loss of energy per unit time be 



\Jlc{x)J 



(A4) 



which has a scaling such as that in GW emission, and the constant 
T is chosen here to be T = 10^^. 

This leads to the following dynamical equations 



Jt+St 



■■Jt+'-\/DjjSt; 



(A5) 
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APPENDIX A: DEPLETION OF PHASE SPACE DUE TO 
GRAVITATIONAL WAVE INSPIRAL 

Throughout this paper, we have not considered the contribution of 
stars in the region where tj > t, (see 92.2. It . Stars in this region do 
contribute to the EMBB: these stars spiral in and become eventually 
directly resolvable EMRIs. Before they can be resolved, however, 
they contribute to a sto chastic background. Thi s contribution to the 
EMBB was studied bv lBarack & Cutler! ( l2004h . and we find that it 
is dominant over the contribution of fly-bys (see figureQ}. 

We stress that our analysis as presented in this paper does 
not apply to the region tj > //, because orbits in this region are 
less likely to be populated. In this appendix we analyse a toy 
model of Monte Carlo (MC) simulations based on the method by 
IShapiro & Mai'chant ( 1978) to sh ow this . For details on the me thod, 
we refer to lShapiro & Marchan3 119781) and lHopmaiil ( |2009|) . The 
model does not include mass-segregation, and is not intended to 
faithfully represent a galactic nucleus, but is instead used to high- 
light the difference in the dynamics in the two regimes delineated 
by tj = tj. 

Consider the following simple dynamical model: 
Energy diffusion — Let the (scaled) energy x perform a ran- 
dom walk such that the diffusion coefficient is 



/)x. = 85x'/4. 



Xt+Sr — Xt 



+ r\/ DxxSt + XGw St, 



(A6) 



where in these equations r is a normally distributed random variable 
with mean zero and unit variance. The three terms in the dynamical 
equations give rise to three time-scales: the time scale for changes 
in j of order of itself due to scattering 



0=fl7// = (l/6.6)x-i/4/ 



(A7) 



the time scale for changes in x of order of itself due to scattering 

t,, = D-Jx^ = ( 1 /85)x" ; (A8) 
and the time for changes in x of order of itself due to GW emission 

(A9) 



_ -'^ _ -4 -7 V/2 
'/ . — X J X^iax'^- 
XGW 



We simulate the dynamical process described above with and 
without GW emission by m eans of a MC simula t ion w ith about 
10^ particles. We follow the IShapiro & Marchan 3 h978l) cloning 
scheme in order to resolve the distribution function over several 
orders of magnitude. We define the distribution f{x,j^) as 



fixj 



2, _ xV4 d^N{x,f) 



rflnxdln p 



(AlO) 



(Al) 



where A' is the normalised number of particles and the prefactor 
x^/4/ p' was added in order to divide out the dependence on energy 
and angular momentum if the distribution is isotropic and there are 
no GWs. We show the resulting normalised steady state distribution 
in figure l lAlb for the case without GW emission and figure l lA2b 
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Figure Al. Distribution f{x,p) of stars without GW emission (arbitrary 
scale). The red soUd line denotes the loss-cone, the blue dotted line delin- 
eates tj =ti, and the green dashed line delineates = (/. 
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Figure A2. Same as previous figure, but now with energy losses to GWs. 
Comparison of the two figures shows that the distribution is depleted 
in presence of GW emission, approximately in the region where (, < 
min(fy,(j) ~ tj. Note that some stars do remain within the tj < mm(tj,tx) 
region; these are inspiraling stars, and th eir contribution to the EMBB was 
accounted for bv lBarack & CufleiUzOoA . 



for the case with GW emission. The figures clearly shows that the 
region where tj <tj is depleted if there are energy losses to GWs. 
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